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We present simultaneous, quality fits to final BRAHMS, PHENIX, 
PHOBOS and STAR data on particle spectra and two-particle Bose-Einstein 
or HBT correlations as measured in Au+Au collisions at y^SNN = 130 
GeV. Using the Buda-Lund hydro model, the best fit is achieved when 
the most central l/8th of the particle emitting volume is superheated to 
T > Tc — 170 MeV. In contrast, we find no indication for such a hot central 
region in a similar analysis of Pb+Pb data at CERN SPS. 

PACS numbers: 24.10.Nz, 25.75.-q, 25.75.Nq, 25.75.Gz 

1. Introduction 

Recently, Fodor and Katz calculated the phase diagram of lattice QCD 
at finite net barion density p. The results, obtained with light quark 
masses about four times heavier than the physical value, indicated that in 
the < /Ub < 200 MeV region the transition from confined to deconfined 
matter is a cross-over, with Tc ~ 172 it 3 MeV which is, within errors, 
independent of the bariochemical potential in this region. At the same 
time, deconfined matter is searched for in Pb+Pb collisions with laboratory 
bombarding energies of 158 AGeV at CERN SPS as well as in Au + Au 
collisions at RHIC, with ^/snn = 200 GeV colliding energies. 



2. Basic properties of the Buda-Lund hydro model 

Here we report on an analysis of the measured single particle spectra 
and two-particle Bose-Einstein correlation functions using the Buda-Lund 
hydrodynamical model [HO]- The key features of this model can be sum- 
marized as follows. 
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1) The particle emitting source is assumed to have a core-halo struc- 
ture ini, where the particle production in the core is governed by local 
thermalization and relativistic, three-dimensional expansion, while particle 
production in the halo is governed by the decays of long-lived hadronic res- 
onances, that have lifetimes larger, or much larger than the characteristic 
5-10 fm/c lifetimes of the hydrodynamically evolving core. In case of pion 
production, such resonances are e.g. uj, r] and r]' . 

2) The freeze-out distribution of the core is described by a three - dimen- 
sionally expanding, locally thermalized source, that reflects the symmetry 
properties of the collision (axial symmetry for central collisions). 

3) The model goes back to known solutions of (relativistic) hydrodynam- 
ics in limiting cases when such solutions are known (e.g. non-relativistic or 
ultrarelativistic limits). 

4) The flow field has a transverse and a longitudinal Hubble fiow, the 
density and the inverse temperature distributions are parameterized keeping 
the means and variances of these distributions only, in such a form that all 
the observables can be calculated analytically from the model. 

5) If the fiow vanishes and the temperature and density are constants, 
the freeze-out temperature can be determined from the single-particle spec- 
tra and the geometrical sizes of the source from two-particle Bose-Einstein 
correlations. However, in the limit of a large, strongly expanding source 
with local temperature inhomogeneities, the two-particle correlations are 
dominated by thermal length scales (sizes within which the Boltzmann fac- 
tor is nearly constant), and the single particle spectra contain (or even are 
dominated) by fiow and temperature inhomogeneity effects, corresponding 
to a scaling limiting behaviour. 

6) In the scaling limit, the analytic formulas given in refs. [21IS1I1] yield 
the analytic prediction that Rout — Rside — Riong ^ ^/y/^t-, or, with other 
words, Rout/ Rside - Rside/ Riong ~ Riong/Rout ^ 1, where the equality is 
reached for sufficiently large transverse masses, satisfying mt ^ Tq, and at 
lower nit the model predicts scaling violating correction terms. 

This Buda-Lund hydro model and the above mentioned analytic formu- 
las were recently reviewed and summarized in ref. as far as the single 
particle spectra and the two-particle Bose-Einstein or HBT radius parame- 
ters are concerned. Since at RHIC the BRAHMS and the PHOBOS collab- 
orations reported data on the pseudo-rapidity distributions in refs. [ZllHl, 
we have extended this model and have improved the analytic calculation in 
such a way, that in the longitudinal, space-time rapidity variable, the posi- 
tion of the point of maximal emittivity can be found exactly now, and thus 
analytic calculation is not anymore limited by the linearized saddle-point 
equations. 
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3. Fits to final run-1 Au+Au data at RHIC 



BudaLund hydro fits to 130 AGeV Au+Au 
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Fig. 1. Solid line shows the simultaneous Buda-Lund fit to final Au+Au data at 
y/sNN = 130 GeV, vers. 1.5. The transverse mass distributions are measured by 
PHENIX PI, the pseudorapidity distributions of charged particles are measured 
by BRAHMS [7| and PHOBOS while the transverse mass dependence of the 
radius parameters are measured by STAR [TI] and PHENIX fTH] , 



BL parameter 


Au+Au RHIC 


Pb+Pb (SPS) 


h+p SPS 


To [MeV] 

Rs\G M 

To [fm/c] 
Ar [fm/c] 
Ar] 

Ts [MeV] 
Te [MeV] 


214 ± 7 
1.0 ± 0.1 
8.6 ± 0.4 
6.0 ± 0.2 
0.3 ± 1.2 
2.3 ± 0.4 
0.5 To fixed 
102 ± 11 


139 ± 6 
0.55 ± 0.06 
7.1 ± 0.2 
5.9 ± 0.6 
1.6 ± 1.5 
2.1 ± 0.4 
131 ± 8 
87 ± 24 


140 ± 3 
0.20 ± 0.07 

7.3 ± 0.8 

1.4 ± 0.1 
1.36 ± 0.02 
1.36 ± 0.02 

82 ±7 
0.1 ± 0.1 



Table 1. Source parameters from simultaneous fits of final BRAHMS, PHENIX, 
PHOBOS and STAR data for Au + Au collisions at ^/sWn = 130 GeV data, as 
given in Fig. 1, obtained with the Buda-Lund hydrodynamical model, version 1.5. 
Pb+Pb data were fitted in ref. T^, while h+p data at CERN SPS in ref. (Tl|. 
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BudaLund hydro fits to 130 AGeV Au+Au 
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Fig. 2. Top row shows the transverse mass dependence of the side, out and 
longitudinal HBT radii, the central line shows their pairwise ratio (usually only 
^out/^side shown) together with the Buda-Lund fits, vers. 1.5. The bottom 
line shows the inverse of the squared radii. The intercept of the curves in this row 
is within errors zero for the two transverse components, so the fugacity is indepen- 
dent of the transverse coordinates. As the intercept is nonzero in the longitudinal 
direction, the fugacity (hence particle ratios) are rapidity dependent. 



In this analysis, final BRAHMS [Z|, PHENIX [31113, PHOBOS and 
STAR fTV data are used, and all these datapoints were fitted simultaneously, 
using the analytic expressions and the CERN Minuit fitting package. The 
fitting package used in this analysis is version 1.5, made public at 6 . A well 
defined minimum was found, with accurate error matrix and a statistically 
acceptable fit quality. The results are shown on Fig. 1 and the essential pa- 
rameters are summarized in Table 1. In addition to the parameters in Table 
1, we have 6 chemical potentials that regulate the absolute normalization 
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of the single particle spectra of vr , K , and p . The quality of the fit to 
the RHIC data is characterized by x^/NDF= 158/180, that corresponds to 
a confidence level of 88 %. More details are given in ref. ^21- Here we note 
only, that the net baryochemical potential, fiB can be determined from the 
difference of the proton and antiproton chemical potentials. In our case, 
this value is hb = 77 it 38 MeV, well within the range where the critical 
temperature is within errors independent of the baryochemical potential 
according to ref. pp. The fitting of RHIC data turned out to be more com- 
phcated than that to the CERN SPS NA44, NA49 and WA98 data, given 
in ref. jl3j . The reason was that the RHIC HBT correlation data are closer 
to the scaling limiting case than the data at CERN SPS. This implied that 
the fugacity distribution, exp{fi(x)/T{x)) becomes a constant parameter at 
RHIC, within resolvable errors, independently of the transverse coordinate 
values. Hence the "geometric" contribution to the HBT radius parameters 
is negligibly small. The finite transverse size seems to be generated by the 
finite transverse size of the hot region, but the fugacity term is within the 
resolvable errors independent of the transverse coordinates in this region. 
This is signalled by the fact that the fits became independent of the geomet- 
rical source sizes, Rq, which characterizes the transverse variations of the 
fugacity term, and we observed that Rq becomes an irrelevant parameter 
at RHIC, with 1/i?^ 0. This can be seen in the lower plots of Fig. 2, 
where we show the inverse squared radii versus m^. The Buda-Lund model 
predicted, see eqs. (53-58) in ref. [2| and also eqs. (26-28) in ^5], that 
the linearity of the inverse radii as a function of rrit can be connected to 
the Hubble flow and the temperature gradients. The slopes are the same 
for side, out and longitudinal radii if the Hubble flow (and the temperature 
inhomogeneities) become direction independent. The intercepts of the lin- 
early extrapolated mt dependent inverse squared radii at = determine 
I/Rq, or the magnitude of corrections from the finite geometrical source 
sizes, that stem from the exp(//(x)/T(x)) terms. We can see on Fig. 3, that 
these corrections within errors vanish at RHIC. This result is important, 
because it explains, why thermal and statistical models can be successful at 
RHIC: if exp(/i(x)/T(x)) = exp(^o/2o), then this factor becomes an overall 
normalization factor, proportional to the particle abundances. Indeed, we 
found that when the finite size in the transverse direction is generated by 
the T(x) distribution, the quality of the fit increased and we had no de- 
generate parameters in the fit any more. In Table 1, we determined Rs at 
RHIC by the condition, that T^r^ = Vy = Rg) = Tq/2. 

On the other hand we have obtained similarly good description of these 
data if we require that the four-velocity field is a fully developed, three- 
dimensional Hubble fiow, with = x^/r, however, we cannot elaborate on 
this point here due to the space limitations |16j. 
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4. Conclusions 

Table 1 and Figures 1 and 2 indicate that there is no real HBT puzzle 
at RHIC- the Buda-Lund hydro model works here as well as at CERN SPS, 
and gives a good quality description of the transverse mass dependence of 
the HBT radii. For the dynamical reason, see refs. and 0. Naturally, 
we find that the transverse size of the source is larger and the transverse flow 
is stronger in heavy ion collisions than in h+p collisions. We also observe, 
that the central temperature, Tq = 214 it 7 MeV, is significantly higher at 
RHIC, than the critical temperature, Tc = 172 it 3 MeV, calculated from 
lattice QCD, and this is a significant, more than 5 standard deviation effect. 
Finding similar parameters from the analysis of the pseudorapidity depen- 
dence of the elliptic flow, it was estimated in ref. [TP that about l/8th of 
the total volume is above the critical temperature in Au+Au collisions at 
RHIC at the time when pions are emitted from the source. We interpret this 
result as an indication for quark deconfinement and a cross-over transition 
in Au+Au collisions at y^s^viv = 130 GeV at RHIC. A similar conclusion 
has been first achieved in ref. ^7] when analysing the final PHENIX and 
STAR data on midrapidity spectra and Bose-Einstein correlations, but only 
at a three standard deviation level. By including the pseudorapidity distri- 
butions of BRAHMS and PHOBOS, the significance of the result increased. 
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